Spin Configuration in the 1/3 Magnetization Plateau of Azurite Determined by NMR 
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High mag netic field 63 ' 65 Cu NMR spectra were used to determine the local spin polarization in the 
1/3 magnetization plateau of azurite, Cu3(C03)2(OH)2, which is a model system for the distorted 
diamond antiferromagnetic spin- 1/2 chain. The spin part of the hyperfine field of the Cu2 (dimer) 
sites is found to be field independent, negative and strongly anisotropic, corresponding to «10 % 
of fully polarized spin in a d-orbital. This is close to the expected configuration of the "quantum" 
plateau, where a singlet state is stabilized on the dimer. However, the observed non-zero spin 
polarization points to some triplet admixture, induced by strong asymmetry of the diamond bonds 
Ji and J3. 

PACS numbers: 75.10.Pq, 75.25. +z, 67.80.dk, 76.60.-k 



The natural mineral azurite, Cu3(C03)2(OH)2, has 
been recently recognized [l| as a model system for the 
frustrated antiferromagnetic Heisenberg spin-1/2 chain of 
"distorted diamond" geometry defined in Fig. [1] Its most 
prominent feature is a large plateau in the magnetization 
curve at 1/3 of the saturation magnetization, which ex- 
tends from 11 to 30 T when the applied magnetic field 
(Hq) is perpendicular to the chains. Such a "1/3 plateau" 
is usually associated with a classical collinear up-up-down 
(uud) type of spin arrangement, or rather to a quantum 
state which has this classical analogue. For example, a 
uud state is predicted for spins on a two-dimensional tri- 
angular lattice and observed in the Cs2CuBr4 compound 
0. The 1/3 plateau in azurite is proposed to be of fun- 
damentally different, "OOu" type, where the dominant 
Ji coupling ensures that the two "dimer" spins on the 
Cu2 sites (see Fig. [T]) are in a singlet state, while the 
third "monomer" (Cul) spin is completely polarized by 
the field. As this state is based on the presence of a 
singlet, it is of pure quantum nature without a classi- 
cal analogue. Azurite is a good candidate to be the first 
system exhibiting such a 1/3 plateau state, but a direct 
experimental evidence is still missing. The point is that 
both types of plateaus are predicted for a diamond chain, 
the OOu type driven by dominant J2 coupling and an ana- 
logue of the uud state in presence of dominant J\ and J3 
[3j,[j|. The azurite is close to the phase boundary between 
them, and there is a controversy on the J values proposed 
from the magnetization, specific heat and neutron scat- 
tering experiments [H, H, [a 0| • The two different plateau 
types are distinguished by very different local spin polar- 
izations, which can in principle be directly accessed by 
performing nuclear magnetic resonance (NMR) on the 
on-site copper 63,65 Cu nuclei. In this letter we present 
such NMR data which show that in the 1 /3 plateau the 
dimer spins are nearly in the singlet configuration and 
thus confirm the 00m type of plateau. We find a small 
non-zero spin polarization of these sites, estimated to ap- 
proximately 10 % of full polarization, which points to an 



important asymmetry of J\ and J3 couplings. The ob- 
served polarization is incompatible with a uud type of 
plateau, in which the dimer spins are strongly polarized. 

In general, the copper NMR spectrum of a single crys- 
tal consists of 6 NMR lines per each non-equivalent Cu 
site, corresponding to three transitions between energy 
levels of a spin 1 — 3/2 nucleus for each of the two 63 Cu 
and 65 Cu isotopes. In the crystallographic structure of 
azurite, shown in Fig. [TJ we recognize two different cop- 
per sites in two equivalent chains of different orientation 
with respect to the arbitrary direction of the applied 
magnetic field Ho (a). We also note that the two Cu2 
sites of each dimer are expected to be undistinguishable 




FIG. 1: (color online) Diamond chains formed by the ex- 
change interactions, JijSi -Sj, between 5 = 1/2 spins of Cu 2+ 
ions in the crystal structure of azurite. There are two equiv- 
alent but differently oriented chains, related by the ac-plane 
of mirror symmetry. Chains contain "monomer" spins on the 
Cul sites coupled by Ji and J3 interactions to each spin of the 
dimer formed by the two Cu2 sites, mutually coupled by J2. 
There is an inversion symmetry on each Cul site and at the 
center of each dimer. For each Cu site 4 nearest neighboring 
oxygen atoms (connected by thin lines) define approximately 
the plane of the local symmetry of the wave functions and of 
the corresponding EFG tensor. Dotted line vectors define the 
rotation angle a, and the angle # between the magnetic field 
and the Z principal axis of the EFG tensor. C and H atoms 
are not shown. 



by NMR, which is ensured by the inversion symmetry 
with respect to the center of each dimer. We thus ex- 
pect that 63 > 65 Cu NMR spectrum has 6 x 2 x 2 = 24 NMR 
lines. Single crystal spectra presented in Fig. [5] contain 
only 12 lines, meaning that we observe only one of the 
two Cu sites. We recall that these spectra are taken in 
the plateau phase, so that they are not affected by the 
Neel ordering appearing at 1.9 K at field values below the 
plateau phase [8J. A standard way for the identification 
of the observed site is to compare the symmetry of the 
local electric field gradient (EFG) tensor determined by 
NMR to what is expected from the local symmetry of the 
four nearest neighboring (NN) oxygen atoms (see Fig.[l}. 
This rather technical procedure, explained in detail in the 
following paragraph, unambiguously demonstrates that 
only Cu2, i.e. the dimer site is observed by NMR. 

Each copper isotope (I = 3/2) generates a triplet of 
NMR lines whose average frequency reflects the Zeeman 
coupling to the total "effective" magnetic field H c ff , while 
the line splitting is induced by the "quadrupolar" cou- 
pling to the local EFG. The corresponding nuclear spin 
Hamiltonian, H = HjI-H eS + /wq [3/| - 1(1+1) + v( I + + 
J— )/2]/6, is uniquely defined by 5 parameters: the EFG 
tensor described by the quadrupolar coupling i/q and 
its asymmetry parameter 77, and H c g and its direction 
(^efg 1 Vefg) with respect to the principal axes (X,Y,Z) 
of the EFG tensor Knowing the gyromagnetic ra- 
tios for the two isotopes, 63 7 and 65 7, as well as the ra- 
tio of their quadrupolar couplings, 63 vq / 65 vq — 1.0805, 
these 5 parameters can be fit to provide the observed 6 
NMR frequencies (for each chain) . The vq and 77 param- 
eters do not depend on the orientation of magnetic field 
and, in particular, they are common to spectra from two 
chains shown in Fig. [1] In this experiment we performed 
in-situ rotation of the crystal around the axis that was 
close to the crystal a-axis, and have taken several com- 
plete spectra at different rotation angles a (see Fig. [TJ, 
as shown in Fig. [2j The smallest NMR line widths and 
thus the most precise fits are obtained when i?efg — 90°, 
which for chain 1 corresponds to a = 45°. This par- 
ticular orientation has therefore been used to determine 
63 j/q = 36.5 MHz and 77 = 0.085, and the obtained EFG 
values have been successfully used to produce all the 
other fits - this time by fitting only three parameters 
(H e g, $efg, </?efg) for each set of 6 NMR frequencies 
(see Fig. [2]). The experimentally obtained i?efg(oO de- 
pendence could then be compared to a simple approxi- 
mate estimate for this quantity, based on the crystallo- 
graphic structure. We know that for an ideal tetragonal 
coordination a pure dx^-Y 2 orbital pointing towards 4 
NN oxygens generates axially symmetric EFG with the 
strongest principal axis along the Z direction. There- 
fore, the best estimate for the Z axis is the normal to 
the plane approximately defined by 4 NN oxygens (see 
Fig. [1]). The direction of H e ff is approximated by the 
direction of the applied field, supposing that the rotation 
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FIG. 2: (color online) Rotation dependence of copper NMR 
spectra of azurite at T — 1.5 K. Magnetic field of 15.0 T is ap- 
plied perpendicular to the rotation axis that was close to the 
crystal a-axis. (a) NMR spectrum taken at u = 85° and the 
corresponding fit (vertical lines), as explained in the text, (b) 
Angular dependence of the observed line positions (crosses, 
with lines to guide the eye) and the corresponding fits (open 
symbols: circles for the central transitions and squares and 
diamonds for the satellites). Color code is given in the figure: 
dark (light) grey lines for the chain 1 (2) and solid (dashed) 
lines for the 63 ^ 65 'Cu isotope. Vertical dotted line denotes 
Ho || ac-plane orientation, where both chains are identical. 



axis is precisely the a-axis of the crystal. These estimates 
of Z and H e a directions define the angle ^Theory, and its 
rotation dependence ^Theory (oc) is plotted in Fig. [3] for 
both Cu sites and both chains, in comparison with the 
experimental $efg(oO values. Neglecting small offset due 
to various approximations, from these data one clearly 
identifies that the observed NMR signal corresponds to 
dimer Cu2 sites and is incompatible with the Cul sites. 

The principal information obtained from the fits is the 
spin part of the hyperfme field, H sp j n = H e ff — (1 + 
K or b) Ho = A gHft(S), induced by the local spin polariza- 
tion (S) through the hyperfme coupling tensor A. The 
orbital (Van Vleck) shift tensor K or b is here a minor 
correction, because typical A values for a copper spin 
are as large as Az ~ — 20T//ie, with large anisotropy 
Az/A± w 10. Knowing that the EFG and the hyper- 
fine shift tensors are dominantly determined by the same 
wave function, we expect that the principal axes of both 
tensors are approximately the same, so that the H sp i n vs. 
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FIG. 3: (color online) Left scale (black lines and symbols): ro- 
tation dependence of the NMR shift obtained from the high- 
est frequency 65 Cu NMR line (open circles), the fit to these 
data (solid line) and an estimate of the orbital shift (dashed 
line) . These results provide an estimate of the hyperfine field 
#ipin ~ —2.0 T and its anisotropy H^ pi JH^, ln » 11, which 
corresponds to 10 % spin polarization in a d-orbital. Right 
scale (color / grey-scale lines and symbols): $efg deduced by 
fitting the line positions shown in Fig. &b) (solid squares and 
diamonds) , following closely the predictions for the dimer Cu2 
site (solid lines) and not the Cul site (short-dashed lines), as 
explained in the text. Experimental $_efg(oO dependence 
(diamonds) is used to define the upper horizontal scale. 



$efg dependence provides a complete information on the 
local spin polarization. In Fig. [3] we have plotted the ro- 
tational dependence of the experimental NMR line shift 
(-ffeff — Hq) and its extrapolation by a sinusoidal fit, to- 
gether with an estimate of the orbital shift (K^ lh w 1.3 %, 
0.3%). From these data we estimate the spin- 
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induced hyperfine field to be -ffj pin 
anisotropy irj[,i n /-ff S p in ~ 11 (where 



=s -2.0 T, with an 
and _L refer to the 



principal axes of this tensor). This corresponds to about 
10% spin polarization of a typical Cu 2+ dx 2 -Y 2 orbital. 
The error in these values is estimated to be ~ 20 %, dom- 
inantly from the extrapolation of the angular dependence 
to H^ pin In particular, in Fig. [3] we clearly see that 
the maximum of the experimental line shift is shifted by 
18° from the expected ^efg = 90° value. This means 
that the principal axes of the EFG and the hyperfine 
tensors are not really parallel, pointing to a departure 
from the simplified picture of pure <ix 2 -Y 2 orbital. In- 
deed, the electronic density observed by x-ray diffraction 
suggest significant admixture of other orbitals 11 1 . 

In a true magnetization plateau the magnetization 
should not vary with the magnetic field. In order to test 
this most prominent feature of the plateau in azurite, we 
performed very high field measurements of the copper 
NMR spectra, in the field range 17-28 T and at 1.4 K, 
for Hq applied close to the c-axis. Because the effects of 
the quadrupolar coupling and the hyperfine shift are en- 



tangled in the NMR spectra, the shift can only be deter- 
mined by the complete NMR fits as explained in the third 
paragraph. The line positions and fits shown in Fig. 0] 
indeed confirm that H sp i n is to a high precision field in- 
dependent in the plateau. For a field variation from 19 
to 26 T (i.e. 37%) the change in the measured |i? sp in| is 
found to be (1 ± 1) %, where the precision is limited by 
our estimate of the orbital shift tensor. This information 
might be important to constrain the possible effects of 
Dzyaloshinski-Moria (DM) interaction terms, which may 
induce some weak field dependence of the spin polariza- 
tions. The DM interaction on the dimer bond has been 
invoked to explain strong anisotropy of the width of the 
plateau [l[. However, the presence of an inversion center 
at the center of the dimer precludes such a term, and 
only DM interaction on J\ and J3 exchange paths are 
possible. Whether or not they can explain the observed 
anisotropy has not yet been studied theoretically. 

Despite considerable efforts, we could not observe the 
NMR signal from the monomer Cul spin. We have tried 
to find it at very low temperature in order to minimize 
the spin fluctuations. In this way the longitudinal (T^ 1 ) 
relaxation rate was reduced, however, the same does 
not necessarily apply to the transverse (T^ 1 ) relaxation 
To provide the correct total polarization of the 1/3 
plateau, the spin polarization at the Cul site has to be 
« 80 %, that is ~8 times more than at the dimer sites. We 
can then roughly estimate that the corresponding T 2 _1 ra- 
tio is of the order of 8 2 = 64. As the dimer T2 values are 
in the 10-100 fxs range (depending on the orientation), 
this factor is enough to reduce the Cul T 2 below the 



3: 



D 

a 

o 

<L> 

tr 




18 19 20 21 22 23 24 25 26 27 28 
Applied magnetic field: H g (T) 



FIG. 4: (color online) Magnetic field dependence of the 
63 ' 65 Cu NMR line positions in azurite at 1.4 K, for the field 
orientation close to the c-axis. The spectra were taken at 
constant frequency by sweeping the field, and results plot- 
ted in the reduced frequency scale to eliminate the dom- 
inant field dependence (7)^0, with (7) = ( 63 7 + 65 7)/2 = 
11.687 MHz/T. Lines are linear fits to the observed field de- 
pendence, and symbols fits to this linear interpolation at 19, 
22.5 and 26 T, which confirm that the spin polarization of the 
dimer Cu2 site is magnetic field independent. Color/symbol 
code is the same as in Fig. [2] 
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experimental dead time for the observation of an NMR 
signal. This provides a reasonable explanation why Cul 
spin could not be observed, but also an important hint on 
the system: the longitudinal spin fluctuations (effective 
in T2 relaxation) are probably not gapped. 

Finally, we remark that here we have only considered 
the standard on-site hyperfine coupling to the copper 
spin, and not the transferred hyperfine coupling which 
could in principle couple the observed nuclear spin at 
the dimer site to the neighboring (strongly polarized) 
monomer spin. This latter mechanism typically relies on 
almost negligible admixture (~1 %) of the on-site s-wave 
orbital in the extended Wannier wave function belonging 
to the neighboring electronic spin. While the on-site spin 
polarization induced in this way is negligible (~1 %), very 
high hyperfine coupling of an s-wave orbital (^200 T//ib) 
can in principle provide significant transferred hyperfine 
field. However, this field is positive and isotropic, in ob- 
vious contradiction to what is observed in azurite. Neg- 
ative and strongly anisotropic iJ sp in necessarily implies 
significant on-site spin polarization of the Cu orbital. 

To provide a simple discussion for the observed po- 
larization of dimer spins, we note that the approxi- 
mate wave functions proposed for the 1/3 plateau of 
the diamond chain Q can be generalized to represent 
an arbitrary mixture of the three single-spin-flip states, 
*(T,/3) = cos/3[ lit?) T|UT)]/\/2- sin/3 |TU>, which 
by construction has correct total and local spin polar- 
ization of and /iBx(sin 2 /3, sin 2 /3, cos 2 /3 — sin 2 /?), 
respectively. In particular, there is equal spin polariza- 
tion on the two dimer spins. In this notation the two 
reference plateau states [3j are 00m = W(— ,0) and uud 
= arccos(l/v3)), where the sign difference corre- 

sponds to the different symmetry. (Note that at least 
two unit cells should be taken into account to properly 
represent all symmetries of the system.) Using ^(t-, 0) a s 
a trial function to minimize the energy by optimizing /?, 
one can easily see that the pure singlet OOu (j3 — 0) state 
is obtained only for the symmetric diamond couplings 
J\ = J3, while deviation from this case necessarily leads 
to some admixture of the triplet, meaning some non-zero 
polarization of the dimer spins (/? ^ 0). Observed po- 
larization on the dimer site means that this admixture 
is significant, sin/? w \/0.1 w 0.3 and thus should 
correspond to an important asymmetry of couplings. A 
correct estimate of the corresponding J coupling values 
should rely on numerical solutions of the spin Hamilto- 
nian, relating the observed spin polarization to the cor- 
responding constraint on the J couplings, say the J3/J1 
vs. J2/J1 dependence [14(. Two available points predict- 
ing the correct local spin polarization, J3/J1 = —0.5 [l3j 
and +0.4 [14|, show that from NMR we cannot directly 
conclude if one of the couplings is ferromagnetic or not. 
However, if NMR data is combined with other constrains, 
as the width of the plateau from the magnetization data 
[H and the energy of excitations from the neutron scat- 



tering data 0, one should clearly define the couplings, or 
indicate whether the diamond chain model is too simple 
to describe azurite. Here we recall a possible influence of 
inter chain couplings. 

In conclusion, by copper NMR in the 1/3 magnetiza- 
tion plateau of azurite we have determined the local spin 
polarization of the dimer spins to be f=a 0.1 /xb- This pro- 
vides the first direct evidence for a "quantum" type of a 
1/3 plateau having no classical analogue, which consists 
of dimers in a singlet state and fully polarized monomers. 
The deviation from ideal zero polarization of the dimer 
implies important asymmetry of the diamond couplings, 
Ji / J3, and provides a strong constraint for the deter- 
mination of their values. 
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